WORLD INTELLECTUAL PROPERTY ORGANIZATION 
International Bureau 




PCT 

INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(51) International Patent Classification ^ : 
GOIR 33A3415, 33/48 



Al 



(11) International Publication Number: WO 99/54746 

(43) InternaUonal Publication Date: 28 October 1999 (28.iO.99) 



(21) International Application Number: PCT/IB99/00652 

(22) International Filing Date: 14 April 1999 (14.04.99) 



(30) Priority Data: 
98201229.6 



17 April 1998 (17.04.98) 



EP 



(71) Applicant {for ail designated States except US): KONIN- 
KLIJKE PHILIPS ELECTRONICS N.V. [NL/NL]; Groe- 
newoudseweg 1, NL-5621 BA Eindhoven (NL). 

(71) Applicant (for SE only): PHILIPS AB [SE/SE]; Kottbygatan 7» 
Kista, S-164 85 Stockholm (SE). 

(72) Inventors; and 

(75) Inventors/Applicants (for US only): PRUESMANN» Klaas, 
P. [DE/NL]; Prof. Holstlaan 6, NI^5656 AA Eindhoven 
(NL). WEIGER. Markus [DE/NL]; Prof. Holstlaan 6, 
NL-5656 AA Eindhoven (NL). SCHEIDEGGER, Markus, 
B. [CH/NL]; Prof. Holstlaan 6, NL-5656 AA Eindhoven 
(NL). BOESIGER, Peter [CH/NL]; Prof. Holstlaan 6, 
NL-5656 AA Eindhoven (NL). 

(74) Agent: SCHOUTEN, Marcus, M.; Prof. Holstlaan 6, NL-5656 
AA Eindhoven (NL). 



(81) Designated States: BR, CN, JP, KR, RU, US, European patent 
(AT, BE, CH, CY, DE, DK, ES, FI, FR, GB, GR, IE, IT, 
LU, MC, NL, PT, SE). 



Published 

With international search report. 

Before the expiration of the time limit for amending the 
claims and to be republished in the event of the receipt of 
amendments. 
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(57) Abstract 

The invention relates to a method for obtaining images by 
means of magnetic resonance (MR) of an object placed in a static 
magnetic field, which method includes simultaneous measurement 
of a number of sets of MR signals by application gradients and an 
array of receiver coils, reconstruction of a number of receiver coil 
images from the sets MR signals measured and reconstruction of a 
final image from a distant dependent sensitivity of the receiver coils 
and the first plurality of receiver coil images. In order to reduce the 
acquisition time the number of phase encoding steps corresponding 
to the phase-encoding gradient is reduced with a reduction factor 
compared to standard Fourier imaging, while a same field of view is 
maintained as in standard Fourier imaging. In this way fast cardiac 
imaging may be possible. According to the invention the calculation 
of complicated matrix inversion can be simplified by determining 
of an image vector of the final image from a combination of a 
generalised inverse of a sensitivity matrix and a receiver coil image 
vector. In this way aliasing artefacts in the final image are reduced. 
Furthermore, the reconstruction method enables application of non 
integer reduction factors. 
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Magnetic resonance imaging method and apparatus. 



PCT/IB99/00652 



The invention relates to a method for obtaining images by means of magnetic 
resonance (MR) of an object placed in a static magnetic field comprising the following steps 
generation of an excitation RF pulse in a pan of the body, 
measurement of a plurality of sets of MR signals using a plurality of receiver 
coils along a trajectory in k-space comprising a first plurality of lines by application of a read 

gradient and other gradients, 

reconstruction of a final image fi-om a combination of the plurality of sets of 
MR signals measured and distance dependent sensitivities of the plurality of receiver coils. 
The invention fiirther relates to an MR apparatus arranged for performing such method. In this 
patent appUcation pixels mean picture elements of a digital image , voxels mean volume-elements 
of a three-dimensional digital object. 



Such a method is known fi-om the article "Simultaneous Acquisition of Spatial 
Harmonic (SMASH): Fast Imaging with radio fi-equency Coil Array" by D.K. Sodickson et aL 
published in Magnetic Resonance in Medicine, vol. 38, page 591-603, 1997. The known 
method is used, for example, in real time cardiac imagmg of human beings. To reduce the 
acquisition time of MR data in the known method a sub-encoding data acquisition scheme is 
used wherein the plurality of sets of MR signals are simultaneously measured using the 
plurality of receiver coils along the trajectory in k-space containing the first plurality of lines 
using the read gradient and the other gradients. The number of lines corresponds to a reduced 
number phase-encoding steps in comparison with a number of phase encoding steps in 
conventional Fourier MR imaging. The receiver coils may be arranged in an array of surface 
coils. A final MR data set is determined fi-om a specific distance sensitivity fiinction of the set 
of receiver coils and the sets of MR signals measured. The final MR data set then contains the 
information of the number of lines of the conventional Founer MR imagmg. The final image 
is then reconstructed by n-ansforming the final MR data set. The number of sets of MR signals 
may be equal to the number of receiver coils in the array. Furthermore, the specific distance 
sensitivity fimction of the set of receiver coils must have a sinusoidal shape. By the 
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simultaneous measurement of MR data is the acquisition time is reduced. The reduction factor 
is determined by the number of lines in k-space corresponding to the final image and the 
number of lines which the sub-encoding data acquisition employs. The reduction in the 
acquisition time may, for example, enable application in real time cardiac imaging or 
functional imaging. A disadvantage of the known method is that it may be limited to only a 
few positions of a field of view and a limited number of orientations of the sHce of the object 
to be image because of the specific sensitivity function of the set of coils. The field of view is 
defined by a distance in a phase encoding direction covered by the trajectory in the k-space. 



It is an object of the present invention to provide an MR imaging method with 
an improved degree of fi-eedom in the choice of field of view and orientation of the slice to be 
imaged. To this end, the method in accordance with the mvention is characterised in that for 
the reconstruction of the final image the method comprises a fiirther step of 

reconstruction of receiver coil images firom each set of MR signals measured 

respectively, and in 

that the final image is reconstructed ft-om a combination of the receiver coil 
images and the distance dependent sensitivities of the receiver coils. In this way the 
reconstruction of the final image does not depend on a specific sensitivity function of the sets 
of receiver coils and can be applied for an arbitrary set of receiver coils and the restrictions on 
the size of the field of view and the orientation of the slice to be image are less severe. 
Funhermore, non-integer reduction factors can be chosen. This option gives a possibility to 
shift aliasing artefacts to less important parts of the image. Simultaneous or partly 
simultaneous measurement of the MR data sets yields a reduction in the acquisition time and a 
possibility of fast MR imaging compared to conventional MR-imaging. 

A particular version of the method according to the invention comprises a step 
of determining an image vector of the final image fi-om a combination of a generalised inverse 
of a sensitivity matrix and a receiver coil image vector, wherein an image vector component 
represents a value of a tissue comrast function at a position of a volume-element selected from 
a first plurality of equidistant volume-elements in a first plurality of adjacent fields of views, 
an element S(id) of the sensitivity matrix represents a sensitivity at the position of the selected 
volume-element with respect to a receiver coil selected from the first plurality of receiver coils 
and a receiver coil image vector component represems a pixel value of a receiver coil image 
corresponding to the selected receiver coil, a position of the pixel in the receiver coil image 
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corresponding with the position of the selected volume-element in the selected field of view. 
The generalised inverse of a matrix S is defined as a matrix product (S"S)-'S", wherem S" 
represents a complex conjugate transpose of the matrix S. The generalised inverse or pseudo 
inverse is known in applied mathematics for the following property 
||5(5*a)-a| =min,||5;t-a|| . 

That is, given a vector equation Sx = a, which does not have an exact solution for x due to 
over determination, the pseudo inverse yields a vector that best fits that equation in the above 
sense. In the reconstruction method according to the invention use is made of this 
minimisation property. The pixel wise reconstruction method provides a possibility to include 
an acmal degree of aliasing in the reconstruction method in order to reduce the different types 
of aliasmg. As a result the method according to the mvention provides substantially aliasing 
free final images m case of both non-integer and imeger reduction factors. 

A fiirther version of the method according to the invention is characterised in 
that a reduction factor of acquisition of the sets of MR signals amounts to a real value smaller 
than or equal to the number of receiver coils. The reduction factor is defined as the factor with 
which the distance between the lines of the trajectory in k-space is increased as compared with 
the distance between adjacent lines in k-space in conventional Fourier MR imaging. As a 
result the acquisition time of the sets of MR signals is reduced proportionally. In cardiac 
imaging, for example, regions of lung tissue of the human being can be often excluded fi-om 
the image reconstruction due to a negligible signal contribution resulting from the choice of an 
optimal reduction factor. An optimal choice of the reduction factor can be chosen such that for 
example a high intensity fold-over of a back chest wall of the human being is directed off the 
heart in to the lung regions. The optimal choice of the reduction factor may result in different 

type of aliasing in the final image. 

A fiirther version of the method according to the invention is characterised in 
that the method comprises a step of determining a dimension of the image vector from a 
number of the first plurality of equidistant voxels in a final field of view of the final image, a 
distance between subsequent voxels being equal to the field of view. An actual degree of 
aliasing can thus be determined and taken into account in the reconstruction process of the 
final image. As a resuh fold-over artefacts in the final image, are reduced. In case of a non- 
imeger reduction factor the degree of aliasmg varies in differem pixels of the final image 
because the acmal number of the volume-elements or voxels varies that contributed to a 
picmre elemem of the single receiver coil images, in general, a pixel value of a single receiver 
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coil image contains contributions from an infinite series of equidistant voxels in adjacent field 
of views. Limiting this series reduces the aliasing in the final image. For example, when the 
reduction factor is 2.5 a first group of pixels of the single receiver coil image will show two- 
fold overiap and a second group of pixels in the single receiver coil image will show tree-fold 
overiap. Contributions of the first group of pixels are determined by two voxels in the final 
field of view, separated by a first distance, and contributions of the second group are 
determined by three equidistant voxels in the final field of view, separated by a second 
distance, wherein the second distance between the voxels is equal to a length of the field of 
view of a single coil image. In the known method either the contributions of two pixels or the 
[) contribution of three pixels can only be taken into account and, as a consequence, three fold- 
over artefacts or twofold fold-over artefact, respectively, are present in the final image. A 
further advantage is that the reconstruction time of the final image is ftirther reduced. 

A ftirther version of the method according to the invention is characterised in 
that the method comprises a step of excluding volume-elements in the final field of view 
5 outside a predetermined region not comprising the object from the reconstruction of the final 
image. In this way a local signal to noise ratio of the final image can be improved. Exclusion 
of a volume-element of a measuremem space from the reconstruction process is possible if 
said volume-element is a priori known to lie outside of the object. In that case no signal 
contribution from this volume-element need to be separated from those contributions of the 
20 volume-elemems of the object in the fold over of the picture elements in the single receiver 
coil images. Knowledge whether a volume-element in the measurement space comprises 
structure of the object can be derived from a first reference image obtained from the object in 
an earlier stage by a conventional MR imaging process, for example an echo-planar imaging 
method. 

25 A ftirther version of the method according to the invemion is characterised in 

that the method comprises a step of determining the sensitivity matrix from a sensitivity map 
comprising a combination of predetermined reference receiver coil images and a 
predetermined reference image. Sensitivity matrices which are independent of tissue contrast 
can thus be obtained for use in the reconstruction method. The obtained sensitivity matrices 
30 reflect only coil sensitivity and allow for smoothing. Said first obtained reference image can 
be used as a predetermined reference image. 

A ftirther version of the method according to the invention is characterised in 
that the method comprises a step of filtering the sensitivity map so as to obtain sensitivity 
maps covering a region in the final field of view larger than the region of interest. A low pass 
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filtered sensitivity map is an example of such a sensitivity map. The imaging method can then 
be applied when the tissue configuration changes, because the sensitivity of the receiver coils 
is in first order independent of the configuration changes with respect to the field of view. This 
enables successive MR imaging of a slightly moving object with a single initial reference 
image measurement, this is advantageous for example, for real time cardiac MR imaging. 
Furthemiore, the sensitivity matrices obtained don't introduce noise in the reconstructed final 
image, so that the signal to noise ratio of the final image is fiuther reduced. A low pass 
filtering step is an example of such a filtering step. Another kind of filtering comprises a step 

of polynomial fitting. 

A fiirther version of the method according to the invention is characterised in 



15 



that the method comprises a step of obtaining the predetermined reference receiver coil image 
by an MR method including the following steps generation of an excitation RF pulse in the 
object, simultaneously measurement of receiver coil reference sets of MR signals by using the 
first plurality of receiver coils along a predetermined trajectory comprising a second plurality 
of Hnes in k-space by application of a read gradient and other gradients and reconstruction of 
the reference receiver coil images firom the receiver coil reference sets of MR signals 
measured, the second plurality of lines being less than the first plurality of lines. A suitable 
knowledge of the coil sensitivities can thus be obtained. For example, a reference 
measuremem of the object is made with the final arrangement of the receiver coils and the 
20 object. A reference image can also be obtained fi-om a combination of the reference receiver 
coil images. 

A fiirther version of the method according to the invention is characterised in 
that the method comprises a step of obtaining the predetermined reference image by an MR 
method compnsing the following sub-steps generation of an excitation RF pulse in the object, 
25 measurement of a reference set of MR signals along a predetermined trajectory comprising a 
third plurality of lines in k-space by application of a read gradient and other gradients and 
reconstruction of the reference image from the measured MR signals measured. This yields the 
predetermined reference image of the object that can be used as the reference image in the 
detemiining of sensitivity matrices. The number of lines of the third plurality may be equal to 
30 the number of lines in k-space corresponding to the final image. The advantage of the 

application of the predetermined reference image thus obtained is that the reconstructed image 
is intensity corrected and contains correct phase information. 

The invention fiirther relates to an MR apparams for obtaining MR comprising 
input terminals for an arrangement of receiver coils, a body coil, means for measurement a 
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plurality of sets of MR signals via the arrangement of receiver coils, means for reconstruction 
of a final image from a distance dependent sensitivity of the receiver coils and the plurality of 
sets of MR signals measured. A fiirther object of the invention is to provide an MR apparatus 
having an improved degree of freedom in the choice of field of view and orienution of the 
shce to be imaged. To this end, the MR-apparatus in accordance with the invention is 
characterised in that the MR apparatus comprises means for reconstruction of a plurality of 
receiver coil images from the plurality of sets of MR signals measured respectively, and means 
for reconstruction of the final image from a combination of the plurality of receiver coil 
images and the distance dependent sensitivities of the receiver coils. 



These and other aspects of the invention are fiirther explained with reference to the 

accompanying drawing of which 

Fig. 1 shows a block diagram of an apparatus for performing a method 

15 according to the invention, 

Fig. 2 shows a block diagram of a detail of such apparatus. 

Fig. 3 shows an example of a final field of view of an object, and 

Fig. 4 shows an example of aliasing in a single receiver coil image. 



20 

An image or distribution may be in two dimensial (x,y) or in three dimensional (x,y,z) 
representation. A digital image consists of picture elements which are usually called pixels or, in 

the three-dimensional case, voxels. 

As is also customary in the art of magnetic resonance imaging, the signals are 
25 sampled and digitised and the fiirther image processing is digital. All that is imphcit in this 
description. 

Where in the description reference is made to nuclear magnetic resonance, this 
imphes nuclear spin resonance which in a specific situation is proton spin resonance and it also 

implies electron spin resonance. 

The object to be imaged is then a patient or generally a human being to be 

examined. 

In Figure 1 a magnetic resonance apparatus 1 is diagrammatically shown. The 
apparatus comprises a set of main magnetic coils 2 for generating a stationary homogeneous main 
magnetic field and several sets of gradient coils 3, 4 and 5 for superimposing on-to the mam 
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magnetic field additional magnetic fields with controllable strength and having a gradient m a 
selecteddirection.ConvenUonally.thedirectionofthemainrnagneticfieldislabelledthez 

direction, the two directionsperpendicular thereto are thexandydirections.Tl.egrad.^^^ 

are energised via a power supply 1 1 • Tl^e apparatus fimher comprises excitation means 6 m the 
form of a so^alled body coil for emittmg radio-frequency pulses (RF-pulses) to an ^ 

7, the excitation means 6 being coupled to modulating means 8 for generating and modulating of 

the RF-pulses. 

Also provided are means tor receiYing the MR-signals, these m«ns can be 
identical to the exciadon means 6 or be separate. a3 with the present invention, m the fonn of an 
„ay of receiver coils 16. A body coil is a coil surrotmding a (par. of a) body and ha. generally a 
very homogeneous sensitivity distnbutton, whereas a receiver coil is a coil used close to one s,d= 
of a (partof a) body. A receiver coil has generally a limited fleld of view, a high sensitivtty but an 
inhomogeneous spatial sensitivity distribution. 

A send-.«:e>ve switch 9 is arranged to separate the received signals fiom the 
excitation pulses. The received magnetic resonance signals are input to teceiving and 
demodulanng means 10. The cxcitatton means 6 and modulating means 8 and the power supply 
1 1 for the gradient coils 3, 4 and 5 arc steered by a control system 12 to generate a predetemtmed 
sequence of RF-pulses and gradient field pulses. The demodulation means ts coupled to data 
processing means 14, for example a computer, for transformatton of the received signals mto an 
1 image that can be made visible, for example on a visual display unit 1 5. 

If the magnetic resonance apparams 1 is put into operation with an object or body 
7 placed in the magnetic field, a small excess of magnetic dipole moments (nuclear or electron 
spins) u, the body will be ah^ed in the direction of the magnetic field, h, equihbdum, this causes 
a net magnedsatton M. m the material of the body 7, dtrected m parallel with the magnenc field. 
5 h, the apparanrs 1 the macroscopic magnedsatton M. ts mampulated by radtattng to the body W- 
pulses having a fiequency equal to dte Larmor Muency of dte dipole moments, thereby bnng,ng 
d,e dipole moments in an exited state and re-orienttng the magnedsation By applymg the 
proper RF-pulses, a rotation of die macroscopic magnetisation is obtained, *e angle of rotation ts 
called the flip-angle. Tire introducdon of variations in the ma^etic field by applying gradten. 
,0 maptetic fields influences dte behaviour of the magnetisadon locally. After the applicatton of RF- 
pulses, the changed magnettsauon will tend to re«m to a state of thetmal equilibrium m the 

magnetic field, eminmg radiatton in Ute process. A well cho» sequence of RF-pulses and 
gradient field puUes causes thtsradtation to be emitted as(decayi„g)ma^edc resonance srgnals 

which ptovide spanal infomtadon about the density and/or relaxation rimes of a^eriam type of 
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nuclei, for example hydrogen nuclei, and the substance in which they occur. By analysis ofthe 
emitted signals, picked up by the coil 6 or with the receiver coil array 16 - and by presentation of 
it in the form of images, information about the internal structure ofthe object or body 7 is 
accessible. The decaying magnetic resonance signal directly after the excitation is called FID 
(Free Induction Decay). By proper application ofthe gradients and/or application of additional 
excitation (RF) pulses also an echo ofthe FID may be produced and used to derive image 
information. 

For a more detailed description of magnetic resonance imaging (MRI) and MRI- 
devices reference is made to the extensive literature on this subject, for example to the book 
"Magnetic Resonance Imaging", by M.T. Vlaardingerbroek, Spinger Verlag, 1996. 

Figure 2 shows the receiving and demodulating means for the array of receiver 
coils are shown. Receiver coils C,,...,q...,C„ arranged in an array of receiver coils 16 are 
connected to receivers circuits R,,...,Rj....Rn analog to digital converters AD,,...,AD,...,ADn, 
memories M,,...,Mj...,Mn, and Fourier transform means FT,,...,FTj...,FTn. By means ofthe 
i Fourier transform means N receiver coil images from MR signals, preferably simultaneously, 
detected by the N receiver coils are determined in a well known way. The N receiver coil 
images are processed in processing means 14 to obtain the final image. In conventional 
Fourier MR imaging MR signals are measured along a trajectory in k-space containing a first 
number NF of parallel lines in k-space by application of appropriate read gradients and other 
20 gradients. In practice this number NF amounts to for example 256. In order to reduce the 

acquisition time ofthe final image a distance between the adjacent lines along which lines the 
k-space have been scanned is enlarged and the number phase encoding steps is reduced 
proportionally. The reduction factor RN can be chosen as a real number or integer number 
smaller than or equal to the number of receiver coils N. In this example the number of phase 
25 encoding steps NF ofthe final image amounts to, for example, 128 and the number of receiver 
coils N amounts to, for example, five. The Fourier transform means reconstructs the single 
receiver coil images from the sets of MR signals, respectively. Each set of MR signal received 
by one ofthe receiver coils. Each single receiver coil image is related to a field of view ofthe 
object which is reduced compared to the final field of view of the object related to the final 
30 image proportionally by the reduction factor. Furthermore, in this example the field of view 
and the final field of view are related to the direction of the applied phase encoding gradients. 

The aliasing artefacts in pixels of a single receiver coil image are caused by a 
superposition of contributions of N equidistant voxels in the final field of view, a distance 
between the subsequent equidistant voxels being determined by the field of view. A fold-over 
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in different pixels a(i) of the receiver coil image vector a having identical positions m the N 
subsequent receiver coil images is described by 

a = Sb 0) 

wherein the vector component a(i) represents a pixel value of a receiver coil image SCI(i) 
selected from the N receiver coil images, an element of the sensitivity matnx S(u) represents 
a complex sensitivity at aposition of the ptxel OPO) selected from the N pixels of the receiver 
coil SC(i) selected from the N receiver coils and the vector b(j) represents a value of a tissue 
10 contrast function at a position of a volume-element in a field of view FOVG) selected from N 
equidistant volume-elements in N adjacent fields of views. The field of view is defined by the 
trajectoiT m the k-space for which the N sets of MR signals of the related receiver coil image 
are measured. The vector b is then determined according to 



15 I Sb-al = rain =>b= S"'a (2), 

wherein S" represents the generalised inverse of the matrix S. ^ ^ ^ 

The generalised inverse of a matrix S is defined as a matrix product (S"S)-'s", wherein S 
represents a complex conjugate transpose of the matrix S. The generalised inverse or pseudo 
20 inverse is knovm in applied mathematics for the following property 
||5(5^a)-al| =min,||Sx-fli 

That IS given a vector equation Sx = a, which does not have an exact solution for . due to over 
detennmation, the pseudo inverse yields a vector that best fits that equation in the above sense. 
In this reconstruction method use is made of this mimmisation property. In equation (2) the 
25 vector b is thus calculated by matnx calculus. These vector b is calculated for all the positions 
.n the field of view. As a result of these steps the final image at the standard resolution can be 

obtained with reduced aliasing artefacts. 

In several cases the aliasing artefacts can be reduced fixrther in the reconstructed 
final image by selecting the reduction factor such that the aliasing artefacts are laying outside 
30 the area of interest. For example, in cardiac imaging region of lung tissue can often be 

excluded from the image reconstruction due a negligible signal contribution. Preferably, an 
optimal choice for the reduction factor will in that case be selected such that a high intensity 
fold-over of the back chest wall is directed of the heart into lung regions. In that case the 



wo 99/54746 PCT/IB99/00652 

10 

optimal choice for the reduction factor may be a non -integer number, smaller than the number 
of the receiver coils N in the receiver coil array 16. However, when the reduction factor 
amounts to a non-integer number different type of aliasing artefacts will occur in the final 
image, for example, when the reduction factor RN has a value 2.5, a two-fold overlap as well 
as three-fold overlap will occur in the final image. 

In order to reduce the different type of aliasing artefacts in the final image an 
aliasing degree can be determined and taken into account by determining a dimension of the 
image vector b. The determining of the actual degree of aliasing and the dimension of the 
image vector b for the receiver coil array 16 comprising, for example three receiver coils, is 
explained with relation to figure 3 and figure 4. 

Figure 3 shows a final field of view of an object. Figure 3 comprises a circle 30 
representing an object to be imaged, for example the heart of a patient or a human being. 
Furthermore, figure 3 shows three adjacent field of views, a cemral field of view 32, a top 
field of view 31 and bottom field of view 33. These adjacent fields of view fit in a final field 
of view 34. If the number of phase-encoding steps of the final image is given by NF and the 
number of phase-encoding steps of the single receiver coil images is given by MR then the 
reduction factor RN is given by RN=NF/NR. Preferably, the reduction factor RN is chosen 
such that an integer number of voxels fits in the field of view 31,32,33 of the single receiver 
coil images. In this example, the number RN equals 2.5. Contributions from voxels of adjacent 
fields of view 3 1,33 and a voxel of a central field of view 32 are then added to yield a pixel 
value of the single receiver coil image. In Figure 3 only the three adjacent field of view 

comprise material of the object 30. 

Figure 4 shows the single receiver coil image 60, comprising a superposition of 
the adjacent field of views 31,32,33. For the reduction of aliasing artefacts in a pixel in a 
receiver coil image of figure 4 the contributing voxels of the final field 34 of view of figure 3 
have to be separated. This superposition of adjacent fields of view is explained by way of 
example. Fig.3 shows three central voxels 35,45,55 in the central field of view 32 in fig. 3 and 
contnbuting voxels 36-39, 46-49,56-59 in the adjacent field of views 31,33. Fig. 4 shows the 
first, second and third pixels 61,62,63 in a single receiver coil image 60, which pixel values 
contains an infinite series of potentially contributing voxels. In fig. 3, for sake of simplicity-, 
only voxels 36-39, 46-49,56-59 are shown. Separation in a single pixel value of the final 
image of the superposition of the pixel values of the receiver coil images is done by use of 
equation (2). In order the reduce the processing time of the final image preferably the 
dimension of the image vector b is determined by the number of the adjacent voxels in the 
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subsequent fields of view 31,32,33, which are located in the final field of view 34, and a 
distance between subsequent equidistant voxels being equal to the field of view. 
Consequently, for a first central voxel 35 close to the edge of the central field of view only two 
voxels 35,36 are present in the final field of view 34 in figure 3, whereas for the second central 
voxel 45 in the centre of the central field of view 32 three voxels 45,46,47 are present in the 
final field of view 34 and for the third central voxel 55 three voxels 55,56,57 are present in the 

final field of view 34. 

In order to improve the local signal to noise ratio of the final image, the method 

comprises a ftirther step of excluding volume-elements in the final field of view outside a 

predetermined region not comprising the object, firom the reconstruction of the final image. 

The signal-to-noise ratio of the reconstructed final image can be approximated by 

?!^=LG4rN (3) 
SNR, 

wherein LG > 1 represents a local geometry factor reflecting a degree of linear dependence of 
the receiver coil sensitivities at superimposed positions. The factor LG is calculated from the 
matrix S only, thus enabling a priori estimates of SRN and optimal choice of RN. In order to 
reduce the geometry factor LG those voxels which He outside of the object as determined by a 
reference image, are excluded from reconstruction and are set to zero in the final image. For 
example, in figure 4 it can be seen that the three voxels 55,56,57 in the final field of view 34 
contribute to the pixel 63 in the single receiver coil image 60. However, from the reference 
image it can be derived that voxel 57 does not represent any structure of the object to 
contribute to the MR signal. So, in this example the dimension of the image vector b 
containing the pixels of the final image is determined by only two contributing voxels 55,56. 
The value of the excluded voxel 57 is set to zero in the final image. 

Different sensitivity matrices S can be determined respectively from sensitivity 
maps of the N receiver coils. In order to obtain the sensitivity maps of the receiver coils a set 
of low resolution reference receiver coil images of arbitrary contrast are divided pixel wise by 
a reference image. A set of low resolution reference receiver coil images are reconstructed, 
) from a set of, preferably, simultaneously acquired MR signals by separate MR imaging 

sequences. The set of simultaneously acquired MR signal can be measured, for example, by a 
known echo planar imaging (EPI) sequence. For example an echo planar imaging sequence as 
described in the cited handbook "Magnetic Resonance Imaging", by M.T. Vlaardingerbroek, 
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Spinger Verlag, 1996 The EPI imaging sequence comprises an excitation of a RF pulse and, 
preferably, simultaneously measurement of five receiver coil reference sets of MR signals 
along a first number of parallel lines in k-space by application of a read gradient and other 
gradients on the static magnetic field. The first number of parallel Unes corresponding with a 
5 reduced number of phase encoding steps compared to the number NF, for example 64 phase 
encoding steps. Preferably, the first numbers of parallel lines are located around a central part 
of k-space such that a low resolution reference receiver coil images are obtained 

The reference image can be obtained, for example, by determining a square root 
of sum of squares of the reconstructed reference receiver coil images. Another possibility to 
1 0 obtain the reference is image is to reconstruct a body coil image from reference MR signals 
measured for example by a separate imaging sequence using the bodycoil 6. Preferably, the 
pulse sequences used to measure the reference MR signals are of the same type as the pulse 
sequences used to measure the receiver coil images. The separate acquisition of MR signal for 
reconstruction of the reference image comprises generation of an RF excitation pulse and 
1 5 measurement of the reference MR signals along a second number of parallel lines in k-space 
by application of read gradient and other gradients, such as a phase encoding gradient. The 
second number is smaller than the number of phase encoding steps NF corresponding to the 
final image, for example, the second number equals 64. Preferably, the second numbers of 
parallel lines are located around a central part of k-space such that a low resolution body coil 
20 image is obtained. Furthermore, the image volume of the reference image may be larger that 
the region of interest of the final field of view. 

In order to enable fast imaging of a slightly moving part of the patient's 
anatomy, such as the patient's heart, with a single initial reference image the method 
comprises a step of filtering the raw sensitivity maps in order to obtain sensitivity maps 
25 covering a region larger than the heart to be imaged. This is advantageous for real time cardiac 
imaging. For example, said step of filtering comprises a low pass filter. In general, sensitivity 
maps may be impaired with noise, in particular in regions with low or zero contribution to the 
MR signals said filter reduces also the noise. 

However, in order to prevent errors at object edges in the sensitivity maps 
30 preferably a two-dimensional polynomial fit technique is applied instead of said low pass 

filtering. According to the two-dimensional fit technique for each pixel of a refined sensitivity 
map a 2D polynomial is locally fit to the sensitivity map. A refined sensitivity value of the 
pixel is obtained by evaluation of the fitting polynomial at the position of the respective pixel 
in the sensitivity map. Fitting is based on valuable data of the sensitivity maps. Valuable data 
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can be obtained by thresholding the sensitivity map according to signal strength in the body 
coil image and subsequent exclusion of isolated pixels, for example, by minimum 
neighbourhood filtering. The remaining pixels, plus a predetermined neighbourhood 
determined by region growing of a region containing the remaining pixels, form the set for 
5 which refined sensitivity values are calculated by fitting. The fitting comprises calculation of 
an L-th order fitting polynomial for the position xo,yo of a sensitivity map of a receiver coil 
given by 



six,y)=Ycy(x-x,ny-yoy 



(4) 

I J=0 



10 



The complex coefficients c-, are determined such that the quantity 
A = Y,^t]s,-s{x,,yti( (5) 

has itl mimmum value, wherein k represents a number of pixels taken into account, denotes 
1 5 the sensitivity value of the pixel xk,yk,taken into account and w^ a weighting factor given by 
the product of a Gaussian kernel centred at XQ.yQ, and the 'reliability' of sk, derived from the 
signal strength in the body coil image applied for the thresholding. A width of the Gaussian 
kernel corresponds to a desired degree of smoothing. The mimmum condition (5) yields 
(l + 1)" linear equations: 



20 



^ = 0 ^ tc,,A,.,=b, (6) 



CCij ,-.y'=o 



wherein 



25 A,,. = l^Mx,-xX'\y.-yor' (7) 

k 



(8) 

Solving the equations (6) for cqq yields the desired refined sensitivity value at the position 

xo.yo- 
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In order to reduce the required calculation time of the polynomial fitting only significant temis 
of the sums in (7),(8) are taken into account and high order fitting is restricted to border 

regions of the objects. 

Furthermore, in order to reduce a scan time the field of view of the final image 
can be chosen so as to be rectangular. For example, in 2D ftmctional brain imaging at equal 
resolution in both directions a 200 x 256 image may be sufficient to measure MR signals with 
EPI imaging sequence comprising 100 phase encoding steps when two receiver coils are 
employed. For a square image of 256x256 more phase encoding steps are necessary with the 
same number of receiver coils. Beside possible applications of the method according to the 
invention in cardiac imaging the described method may also be applied in for example, 
functional imaging of the brain of a human or animal body. 
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CLAIMS: 



1 . A method for obtaining images by means of magnetic resonance (MR) of an 
object placed in a static magnetic field comprising the following steps 

generation of an excitation RF pulse in a part of the body, 
measurement of a plurality of sets of MR signals using a plurality of receiver 
5 coils along a trajectory in k-space comprising a first plurality of lines by application of a read 

gradient and other gradients, 

reconstruction of a final image firom a combination of the plurality of sets of 
MR signals measured and distance dependent sensitivities of the plurality of receiver coils, 
characterised in that for the reconstruction of the final image the method comprises a further 
10 step of 

reconstruction of receiver coil images firom each set of MR signals measured 
respectively, and 

that the final image is reconstructed firom a combination of the receiver coil 
images and the distance dependent sensitivities of the receiver coils. 

15 

2. A method according to Claim I characterised in that for reconstruction of the 
final image the method comprises a step of determining an image vector of the final image 
fi-om a combination of a generalised inverse of a sensitivity matrix and a receiver coil image 
vector, wherein 

20 an image vector component represents a value of a tissue contrast fimction at a position of a 
volume-element in a field of view selected fi-om a first plurality of equidistant volume- 
elements in a first plurality of adjacent fields of views, 

an element S(ij) of the sensitivity matrix represents a sensitivity at the position of the selected 
volume-element with respect to a receiver coil selected from the plurality of receiver coils and 
25 a receiver coil image vector component representing a pixel value of a receiver coil image 
corresponding to the selected receiver coil, a position of the pixel in the receiver coil image 
corresponding with the position of the selected volume-element in the selected field of view. 
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3. A method as claimed in Claim 1 or 2, wherein a reduction factor of acquisition 
of the sets of MR signals amounts to a real value smaller than the number of receiver coils. 

4. A method as claimed in Claim 1 , 2 or 3, wherein the method comprises a step 
of determining a dimension of the image vector from a number of equidistant voxels in a final 
field of view of the final image, the number being smaller than or equal to the number of 
receiver coils, a distance between subsequent voxels being equal to the field of view. 

5. A method as claimed in Claim 1, 2, 3 or 4, wherein the method comprises a step 
of excluding volume-elements in the final field of view outside a predetermined region not 
comprising the object from the reconstruction of the final image. 

6. A method as claimed in Claim 1, 2, 3, 4 or 5, wherein the method comprises a 
step of determining the sensitivity matrix from a sensitivity map comprising a combination of 
a predetermined reference receiver coil images and a predetermined reference image. 

7. A method as claimed in Claim 5 or 6, wherein the method comprises a step of 
filtering of the sensitivity map so as to obtain sensitivity maps covering a region in the final 
field of view larger than the region of interest. 

8. A method as claimed in Claim 6 wherein the method comprises a step of 
obtaining the predetermined reference receiver coil image by an MR method including the 
following steps 

generation of an excitation RF pulse in a part of the body, 

measurement of a plurality of receiver coil reference sets of MR signals using 

the pluraHty of receiver coils along a predetermined trajectory comprising a second plurahty of 

lines in k-space by application of a read gradient and other gradients, and 

reconstruction of reference receiver coil images from the plurality of sets MR 

signals measured respectively, the second plurality of lines being less than the first plurality of 

lines. 



9. A method as claimed in Claim 6 or 7 wherein the method comprises a step of 

obtaining the predetermined reference image by an MR method comprising sub-steps of 
generation of an excitation RP pulse in the object, 
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measurement of a reference set of MR signals along a predetermined trajectory 
comprising a third plurality of lines in k-space by application of a read gradient and other 
gradients and 

reconstruction of the reference image from the measured MR signals. 

10. An MR apparatus for obtaining MR images of a part of a body using a method 

as claimed in anyone of the claims 1-9 comprising input terminals for an arrangement of 
receiver coils, a body coil, means for measurement of a plurality of sets of MR signals via the 
arrangement of receiver coils, means for reconstruction of a final image from a distance 
dependent sensitivity of the receiver coils and the plurality of sets of MR signals measured, 
characterised in that the MR apparatus comprises means for reconstruction of a plurality of 
receiver coil images from the plurality of sets of MR signals measured respectively and means 
for reconstruction of the final image from a combination of the plurality of receiver coil 
images and the distance dependent sensitivities of the receiver coils. 
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